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SURVEY  OF  EUROPEAN  STUDIES  ON  THE  FAR  FIELD 
CHARACTERISTICS  OF  A  SOURCE  OSCILLATING  AND 
TRANSLATING  NEAR  A  FREE  SURFACE 


1.  INTRODUCTION 

Starting  from  the  late  1940’s  and  continuing  to  the  present  day,  a  number  of  significant  papers 
have  been  written  in  the  European  languages  (French,  German,  and  Russian)  investigating  various 
aspects  of  the  Green’s  function  for  a  pulsating  three-dimensional  source  which  is  translating  in  a 
straight  line  below  the  free  surface  of  an  infinite  fluid.  This  Green’s  function  constitutes  the  basic 
building  block  for  many  modern  day  computer  programs  which  analyze  the  motions,  forces,  and  wave 
patterns  of  ships  in  a  seaway.  In  addition,  this  Green’s  function  contains  as  two  special  limiting  cases 
the  widely  used  functions  for  zero  frequency  (Kelvin  wake  case)  and  zero  forward  speed  (circular  ring 
wave  case). 

This  report  summarizes  the  leading  papers  written  in  French,  German,  and  Russian  on  the  for¬ 
mulation,  construction,  analysis  of  wave  pattern,  and  approximate  evaluation  of  wave  amplitude  of 
this  important  Green’s  function.  Since  many  of  the  papers  use  similar  methods  in  treating  the  above 
topics,  the  approach  taken  here  is  not  to  present  a  detailed  summary  of  each  paper.  Instead,  for  each 
topic,  the  report  describes  in  some  detail  the  major  approaches  used  in  these  papers. 

On  the  whole,  the  notation  used  in  the  equations  follows  that  of  the  original  papers.  However, 
for  the  sake  of  uniformity,  the  original  equations  have  been  rewritten  (when  necessary)  to  conform  to 
the  coordinate  system  shown  in  Fig.  1  and  the  consistent  notation  for  principal  variables,  shown 
below.  Figure  1  shows  that  the  x  axis  is  in  the  direction  of  motion  of  the  source.  In  those  studies 
where  the  source  is  taken  to  be  fixed,  the  x  axis  then  coincides  with  the  direction  of  the  current.  The 
y  axis  is  the  horizontal  axis  perpendicular  to  x ,  so  that  the  xy  plane  represents  the  undisturbed  free 
surface.  The  z  axis  is  directed  vertically  upwards. 

Principal  Notation 

k  wavenumber 

r  =  ^{x  -  x')2  +  (y  -  y')2  +  (z  ~  z')2 

r,,  r'  =  n/(jT^-T')2  +  (y  -  y')2  +  (z  +  z')2 

R  =  ^(x  -  x')2  +  (y  -  y')2 

U  speed  of  source  or  current 

x  ,y  ,z  field  point  location 

x',y',z'  source  location 
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(3,y  =  tan  1  0  -  y')/(x  -  x '),  spatial  direction 

6,u  wave  propagation  direction 

v  =  Uu/g ,  ratio  of  source  speed  to  phase  speed 

03  circular  frequency 

2.  CONSTRUCTION  OF  GREEN’S  FUNCTION 
2.1  Basic  Formulation 

It  is  well  known  that  the  Green’s  function  G(x,y,z,t)  for  an  oscillating  trans 
located  at  (x',y',z')  should  satisfy  the  following  conditions  in  a  fixed  coordinate  system 

V2G  =  0 ;z  <0;  (x ,y,z)  *  (x’,y\z') 


d2G  ,  „  ,  dG(x,  y,  z,  t) 

—y(x,y,z,t)  +  g  - -  =0  on  z  =  0 

dt 1  oz 


V  dG  o 

lim  —  =  0 

Z - Oo  OZ 


lim  VG  =  0.  (2. Id) 

«-0= 

In  addition,  a  radiation  condition  of  no  incoming  waves  at  R  —  oo  is  needed  to  ensure  uniqueness  of 
the  solution. 

In  a  coordinate  system  moving  with  the  velocity  U  in  the  x  direction,  the  only  explicit  change  is 
that  Eq.  (2.1b)  is  replaced  by 

d2G  a2c  .  rr7  a2G  dG  n 

■^J-wW+u  +*  fc7*0on:,  *°  ,22' 

where  the  subscript  1  denotes  independent  variables  for  a  moving  coordinate  system.  This  follows 
from  the  well  known  kinematic  relationship  between  time  derivatives  in  fixed  and  moving  coordinate 
systems 


source 

a 

(2.1a) 

$ 

& 

(2.1b) 

i 

(2.1c) 

s#* 

vS-' 

(2. Id) 

v>5 

a  a 

dt  dt  | 


dt 2  dt} 


4  -  2u  ♦  v‘  4 

3/ 12  ar.ajc,  a.tt 


(2.3a) 


(2.3b) 


Grekas  [2J  shows  that  the  free  surface  equation  (2.2)  also  applies  in  a  fixed  coordinate  system  for  a 
fixed  source  in  the  presence  of  a  uniform  stream  U  in  the  -x  direction. 
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The  construction  of  the  Green’s  function  in  most  studies  uses  the  Fourier  transform  approach. 
Here,  the  description  largely  follows  the  lucid  and  detailed  approach  given  by  Bougis  [1].  In  particu¬ 
lar,  he  gives  the  key  differences  between  approaches  using  moving  and  stationary  coordinate  systems. 
Others  who  have  used  similar  approaches  include  Grekas  [2,3],  Haskind  [4],  and  Sretenskii  [5].  The 
approach  by  Sretenskii,  which  gives  a  Fourier  transform  for  the  wave  elevation  itself,  is  also  briefly 
described.  The  convolution  time  domain  approach  used  by  Brard  [6,7]  is  described  in  the  next  sec¬ 
tion.  The  approach  of  Simmgen  [8],  which  obtains  G  as  a  distribution  of  dipoles  and  is  given  in 
single  integral  form,  is  described  in  the  last  section. 

2.2  Fourier  Transform  Approach 

As  is  usual  in  this  type  of  approach,  Bougis  expresses  G  in  the  fixed  coordinate  system  as  the 
sum  of  two  terms 

G  -  —  +  Gx(x,y,zs',y',z',t)  (2.4) 

where  the  first  term  is  a  Rankine  source.  Noting  that  the  elemental  wave  solution 

W(x,y  ,z)  =  cos  9+y  sin  e)  (2.5) 

satisfies  Eqs.  (2.1a)  and  (2.1c),  G  is  written  in  the  following  Fourier  transform  form  over 
wavenumber  k  and  wave  direction  6  space 

1  I  t/2  oo  I 

Gj  =  —  Re  j  j  2  dd  JQ  Gx(6,k ■,t)elaeik(x  cos9+y sin  k  (2.6) 

For  the  sake  of  brevity,  the  use  of  a  fictitious  viscosity  coefficient  by  Bougis,  which  serves  to  deter¬ 
mine  the  proper  sense  of  the  contour  of  integration  around  the  poles  to  satisfy  the  radiation  condition 
but  which  is  ultimately  set  equal  to  zero,  is  omitted  in  the  present  description.  The  use  of  the  ficti¬ 
tious  viscosity  is  described  in  greater  detail  in  Chapter  4,  in  connection  with  Haskind’s  work  [4], 
Using  the  key  identity 


1  =  1  Re 
r  tc 


X/ Z  00 

(  de  f  e~kil  ~l">eikwdk 

J-ir/Z  J0 


where 


w  =  (x  -  x')  cos  6  +  (y  -  >•')  sin 
and  noting  that  in  a  fixed  coordinate  system 


(2.7a) 


(2.7b) 


the  free  surface  condition  (2.1b)  leads  to  the  following  differential  equation  for  the  transform 

Gx(e,kj) 

"a2  "1 

— —  +  gk  Gx(d,k, l)  =  (o>2  +  U2k2  cos2  6  +  gk)  cos  - 

Jt  J 
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(2.9) 


2iUk  cos  6  03  sin  o>t 


,kz' ^ik{x'  cos$+y'  sin  9)  J_ 

k 


Bougis  points  out  that  it  is  not  convenient  to  solve  this  equation  since  x'  is  a  function  of  time.  By 
expressing  cos  oit  and  sin  o3t  in  terms  of  complex  exponentials  and  the  following  transformation 
between  moving  (xj.y^Zi)  and  fixed  ( x,y ,z )  coordinate  systems 


x  =  x,  +  Ut 


(2.10a) 


y  =  Vi 


(2.10b) 


z  =  Z| 


(2.10c) 


Eq.  (2.9)  becomes 


JL  +  aJk 

at2  8 


A  /a  ,  .  1  kz',  -ik(x‘.  cos  ®  +  v;  sin  0) 

G  ,(0,*;O  =  —  e  e 


ei(o,-Vk  cos  9)i  [w2  +  V2k2  cqs2  0  +  gk  -  2Vk  COS  6  03]  + 


e  ,(w+£/*  cos9)'  [cj2  +  U2k2  cos2  0  +  gfc  +  2t/fc  cos  6  w]  k 


(2.11) 


Solving  this  equation  for  G,(0,fc;f),  discarding  the  transient  solutions,  again  using  the  key  iden¬ 
tify  (2.7),  and  returning  to  the  fixed  coordinate  system  by  using  Eqs.  (2.10),  the  following  form  for 
G  is  obtained 


c  =  E 


±  COS  03t 


(x  -  x'Y  +  (y  -  y')2  +  (z  =f  z')i 


.'x2,\/2 


1 


Re 

IT 


*/2 


e*'"'  f  f 

j-t/2  j0 


.*(:  +:'  +m) 


(w  =F  l/fc  cos  0)~ 

where  the  terms  are  summed  for  the  upper  and  lov/er  signs. 


dk 


(2.12) 


Sretenskii  [51  uses  a  similar  approach  to  obtain  G  and  proceeds  further  to  obtain  an  expression 
for  the  wave  elevation  f  in  Fourier  integral  form.  He  considers  a  fixed  coordinate  system  with  flow 
V  in  the  +x  direction  past  a  fixed  pulsating  source  located  at  x'  =  0,  v'  -  0,  with  depth  of  submer¬ 
gence  z'.  Here,  the  sign  of  the  term  involving  2 U  in  the  free  surface  condition  (2.2)  is  changed  from 
minus  to  plus,  resulting  in 


d2G  dG 


dr 


dxdt 


Ox 


,  +  g  -r-  =  0  °n  z  =  0. 
L  dz 


(2.13) 


As  in  the  preceding  formulation,  G  is  taken  to  be  the  sum  of  a  Rankine  source,  its  mirror  image 
about  the  free  surface,  and  a  harmonic  function,  as  follows 


#.•  h'.ii*  ***  ♦,*  i> 


G  =  —  — - —  cos  o>t  +  G'  cos  a )t  +  G"  sin  a >t 

4x  r  r 


(2.14) 


where  0  is  the  strength  of  the  pulsating  source.  By  defining  G]  as  the  complex  sum  of  G'  and  G" 


G,  =  G'  +  iG" 


(2.15) 


substituting  Eq.  (2.14)  into  the  free  surface  condition  (2.13),  and  using  the  key  identity  (2.7a),  G,  is 
obtained  in  the  following  double  integral  form  similar  to  the  last  term  of  Eq.  (2.12) 


G,  =  (  °°  kek(l~l)dk  C* 

1  4x3  Jo 

=  f”  kek(z~l  )dk  f'  - 

a-3  Jq  p 


gik( x  cos  B+y  sin  6) 
gk  ~  (u  +  Uk  COS  d)2 


,ikR  cos  (B-y) 


4t3  j°  gk  ~  («  +  Uk  COS  0)2 

where  use  has  been  made  of  the  polar  coordinate  relations 

x  =  R  cos  7,  y  =  R  sin  y. 

The  wave  elevation  f  is  given  by 


(2.16) 


(2.17) 


f  *  Re  io>G,  +  U 


e,ut  at  z  =0. 


(2.18) 


Substituting  Gt  given  by  the  second  equality  in  Eq.  (2.16)  gives  f  in  the  following  double  integral 
form 


.  _  Re  f  igQ  rT  r“  k{ u  +  Uk  cos  fl)  e  faV**  cos(8  y) 
I  4-r2  o  gk  -  (u  +  Uk  cos  0)2 


dk  r  e 


(2.19) 


2.3  Convolution  Time  Domain  Approach 


Brard  [6]  starts  with  the  two-dimensional  case  of  a  doublet  at  (x'  ,z')  moving  in  an  arbitrary  path 
parallel  to  the  vertical  xz  plane  and  with  strength  M  varying  as  a  function  of  t .  Letting  v  =  x  +  iz, 
v'  =  x'  +  iz',  the  potential  of  a  doublet  in  an  infinite  fluid  is 


FQ{v,t)  =  “  -Qp  l°g  (v  -  v')Af(r)  =  e'*(v  y  )M(t)dk. 


Expressing  G  as  the  sum  of  F0  and  F 


G  =  Fn  +  F 


and  further,  F  in  terms  of  its  Fourier  transform  P 

F(v,f)  =  j"  e  ~ik',P(k  ,r  )dk 

the  free  surface  condition  (2.1b)  leads  to  the  following  differential  equation  for  P 


(2.20) 


(2.21) 


(2.22) 


H 

n 

K 


SR 

m 


Kg 

IS 

m 


fd2  ' 

r  ~\ 

d 2 

~rj  +  8k 
dr 

K.  J 

P  =  - 

v.  J 

e‘kv  M 


(2.23) 


where  a  bar  denotes  the  complex  conjugate.  The  solution  of  this  equation  for  P  and  its  substitution 
into  Eq.  (2.22)  results  in  the  following  final  expression  for  F  (\,t)  in  double  integral  form 


F(v,r)  =  \°°  eik(v -v,)M(t)dk 
jo 

-2  j“  dk  e-'*v  j'  j-  (eikV'M)  cos  ^gk  (t  -  r)dr  (2.24) 

where  the  first  term  on  the  right  hand  side  corresponds  to  an  image  doublet  which  is  symmetrically 
located  above  the  free  surface. 


For  the  generalization  to  the  three-dimensional  fixed  coordinate  system  (x.v.z).  the  doublet  term 
F0  is  replaced  by  the  source  with  strength  S  ( t ) 

•Ao  =  — ^  ek{z-')J0(kR)S(t)dk 


=  ~  j"  dkS(t)  eik[w+iiz+:  )]de  (2.25) 

where  w  is  defined  in  Eq.  (2.7b).  Comparing  Eqs.  (2.20)  and  (2.25)  shows  that  the  generalization 
from  two  to  three  dimensions  involves  an  additional  integration  of  the  wave  propagation  direction  9 
from  0  to  2x,  and  the  following  two  related  substitutions: 

v  =  x  +  iz  x  cos  0  +  y  sin  9  +  iz  (2.26a) 


v'  =  x’  +  iz' 


x '  cos  6  +  v '  sin  9  +  iz 


(2.26b) 


The  resulting  expression  for  <i>,  the  three-dimensional  equivalent  of  F,  the  regular  part  of  the  Green’s 
function,  then  becomes 


2t 


4>(x,y,z,r)=  +  —  f°°  dk  f  d9e~* 

r  7T  J0  Jo 


ik(x  cos  8  *  \  sin  8  +  iz ) 


('  A 

Jo  dr 


gik |ji’(t)  cos  0  +  y'(r)  sin 


|cos  [v'g*  (t  -t)]  Jdr  (2.27) 

where,  analogous  to  the  two-dimensional  case,  the  first  term  represents  an  image  source  symmetri¬ 
cally  placed  above  the  free  surface. 

The  above  expression  represents  the  Green’s  function  for  a  source  with  strength  SU)  and  loca¬ 
tion  x'(t),  y'(t),  z'(t)  varying  arbitrarily  with  time.  To  specialize  to  the  case  considered  m  the 
present  report  of  a  sinusoidally  pulsating  source  moving  rectilinearly,  such  that 


6 


AjVIl 


».* 


I«,ll, 


$  becomes 


jc '(f)  =  f/r,  3; ' (r )  =  0,  z'(t )  =  const  <  0,  S(f)  =  5  cos  wt 


*(x,y,Z,t)  =  —  +  —  \2V  dk  [2r  e~ik[x  cos  #+>  sin  9+i(z+z'nde  x 

r'  4t  Jo  Jo 


(2.28) 


(  i(kU  cos  0  +  w)  ei(kUcos  9+u)t  +  i  (ikC/cos  0  -  w)  <?,(t{/cos  ®-")T 

•-QD 


g-i'fgk.T  gi'Tgkt  +  ei'fgkr  g-i'fgk! 


d  t. 


(2.29) 


2.4  Superposition  of  Dipoles 

Simmgen  [8]  carries  out  his  formulation  in  a  fixed  (x,y,z)  coordinate  system,  with  a  uniform 
flow  U  in  the  jc  direction.  The  double  integral  key  identity  in  Eq.  (2.7a)  is  rewritten  in  single 
integral  form,  as  follows 


where 


l 

r 


J_ 

T 


(2.30) 


w  =  (x  -  x’)  cos  a  +  (y  -  y')  sin  a  +  i(z  -  z’)  =  p  +  i(z  -  z’).  (2.31) 

Since  Re  (i /w)  represents  a  dipole  in  the  complex  pz  plane,  Eq.  (2.30)  shows  that  a  Rankine  source 
may  be  represented  by  the  integral  over  the  inclination  a  of  a  distribution  of  dipoles.  Similar  to  the 
case  of  a  source,  the  presence  of  a  free  surface  for  a  pulsating  dipole  ga  requires  the  addition  of  a 
harmonic  function  ka,  as  follows 


ga  =  Re  —  +  ka  (2.32) 

w 

where  the  time-dependent  multiplicative  factor  e"*  has  been  omitted.  The  Green’s  function  G2  is 
then  obtained  as  the  integral  over  a  of  ga 

C2  =  -  !'*  (2.33) 


The  harmonic  function  ka  is  determined  by  the  free  surface  condition  in  the  form 
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ko  3 

- - —  v  cos  a.  +  — 

cosza  op 


lo 


cos  a 


dz 


ga  =  0  on  z  =0 


(2.34) 


where  v  =  Uu/g  and  k0  =  g/U2.  Adopting  an  approach  used  by  Haskind  [10]  to  study  the  Green’s 
function  for  the  two-dimensional  case  of  a  pulsating,  translating  source,  ga  is  expressed  in  terms  of 
the  functions  F  * 
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1  r  . 

1 
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W, 


8a  =  Re 
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where 
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where 


\(y)  =  v  +  —  cosh  7  +  p(^)  cosh  7 

p<7>  =  cosh27  +  v  cosh  7  . 

The  paths  Lx  and  L2  lie  in  two  sheets  of  the  7-plane,  as  shown  in  Fig.  2,  which  is  adapted  from  Fig. 
3  of  Simmgen.  Eggers  [9]  points  out  that  70  in  this  figure  is  missing  the  term  i  tt/2 

7o  =  x  +  it  +  i  t/2  .  (2.42) 

Otherwise,  the  symmetry  between  +(y  —y')  and  —  (y  —y')  would  be  lost. 

3,  FAR  FIELD  WAVE  PATTERNS 

Bougis  [1],  Grekas  [2,3],  Haskind  [4],  Sretenskii  [5],  Brard  [6,7],  Becker  [11],  and  Eggers  [12] 
study  the  far  field  wave  patterns  due  to  the  Green’s  functions  derived  in  the  preceding  chapter.  With 
the  exception  of  [2,3,4],  which  are  discussed  in  the  following  chapter,  the  investigations  are  limited  to 
obtaining  the  far  field  curves  of  constant  phase,  for  example,  the  wave  crests.  Becker  [11]  points  out 
that  such  a  presentation  may  be  considered  to  be  an  instantaneous  photograph  of  the  free  surface  with 
neighboring  wave  crests  differing  by  a  phase  of  2ir. 

There  are  two  principal  methods  for  determining  the  curves  of  constant  phase.  The  most  com¬ 
mon  approach  is  to  use  the  well  known  method  of  stationary  phase  [1,2, 3, 5, 6,7].  For  the  Green’s 
functions  obtained  by  using  the  Fourier  transform  approach  the  description  here  will  largely  follow 
the  concise  and  direct  approach  used  by  Sretenskii  [5J.  His  unfortunate  error  in  choosing  the  integra¬ 
tion  paths,  which  leads  to  an  incorrect  placement  of  his  wave  patterns,  is  clearly  pointed  out.  The 
wave  patterns  shown  are  those  recently  obtained  by  Bougis  [1],  A  description  is  also  given  of  the 
application  of  this  method  by  Brard  [6,7]  to  his  convolution  time  domain  Green’s  function,  Eq. 
(2.29).  The  other  method  is  the  approach  taken  in  [11,12],  where  the  wave  patterns  are  constructed 
from  basic  considerations  concerning  the  frequencies  of  the  waves  observed  in  moving  and  stationary 
coordinate  systems.  This  may  be  considered  a  more  physical  alternate  approach  as  compared  to  the 
more  mathematical  method  of  stationary  phase.  The  description  here  largely  follows  the  formulation 
given  by  Becker  [11].  The  principal  differences  between  the  approach  taken  by  Eggers  [12]  from  that 
of  Becker  [11]  are  also  pointed  out. 

3.1  Method  of  Stationary  Phase  Applied  to  Fourier  Transform  Approach 

The  usual  procedure  is  to  consider  the  double  integral  in  the  various  expressions  for  Green's 
functions  (Eqs.  (2.12),  (2.16),  or  (2.19))  as  a  single  integral  of  the  form 

/  =  J  h(0)e,Rfmdd. 

For  large  R ,  the  method  of  stationary  phase  states  that  /  may  be  approximated  by 


(3.1) 


where  60  is  given  by 


f'(0o)  =  0.  (3.3) 

If  there  is  more  than  one  root,  i.e.,  point  of  stationary  phase,  then  Eq.  (3.2)  must  be  summed  over  all 
these  roots.  The  entire  Eq.  (3.2)  must  be  considered  if  one  is  interested  in  the  wave  amplitudes. 
This  would  involve  an  evaluation  of  the  integral  over  k  in  order  to  obtain  h(80).  However,  if  only 
the  curves  of  constant  phase  are  of  interest,  then  only  Eq.  (3.3)  need  be  considered.  From  either  Eq. 
(2.16)  for  the  Green’s  function  G\  or  Eq.  (2.19)  for  the  wave  elevation  f,  f  (8)  is  given  by 

f  (6)  =  k  cos  (6  -  7)  (3.4) 

where  kl2  are  the  roots  of  the  denominator  of  the  integrand  in  Eqs.  (2.16)  or  (2. 19) 


gk  -  (a  +  Uk  cos  0)2  =  0 


(3.5a) 


and  are  given  by 


*1.  k2 


(g  —  2  Ua  cos  8)  ±  ^g(g  —  4 Ua  cos  8) 
2  U2  cos2  8 


(3.5b) 


Eqs.  (3.5)  arise  due  to  the  fact  that  the  residues  of  the  inner  integral  in  Eqs.  (2.16)  or  (2.19)  are 
evaluated  at  the  poles  given  by  Eq.  (3.5b).  Taking  the  derivative  of  f  (6)  in  Eq.  (3.4)  and  setting  it 
equal  to  zero  results  in 


sin  (8  -  7)  _  -21) 
sin  7  g 


(a  +  Uk  cos  8). 


(3.6) 


Substituting  the  values  of  k,,fc2  into  Eq.  (3.6)  gives  two  equations  relating  the  wave  propagation 
direction  8  to  the  spatial  direction  7 


^  _ 

-ctn  7  =  tan  6  ±  - —  ^  -  Av  cos  8  (3.7) 

sm  2d 

where  the  upper  sign  is  for  k{,  the  lower  sign  for  k2 ,  and  v  =  Uu/g .  The  term  -4  v  cos  8  under 
the  square  root  sign  suggests  that  the  characteristics  of  the  wave  patterns  may  be  significantly  dif¬ 
ferent  for  v  ^1/4. 

Since  Eq.  (3.7)  holds  for  either  G[  or  f,  this  means  that  the  far  field  curves  of  constant  phase 
are  identical  for  both  G \  and  f .  The  actual  amplitudes  are,  of  course,  different  since  the  integrands 
are  different  in  Eqs.  (2.16)  and  (2.19). 

The  proper  ranges  of  8  to  be  used  in  Eq.  (3.7)  are  determined  by  placing  the  contours  of 
integration  in  the  inner  integral  with  respect  to  k  above  or  below  the  poles  k ,  and  k:  to  satisfy  the 
radiation  condition  that  the  energy  of  all  the  wave  systems  move  away  from  the  source  [l].  Sretenskii 
chooses  paths  of  integration  which  surround  all  the  poles  with  upper  semicircles.  Unfortunately,  this 
leads  to  systems  of  divergent  and  transverse  waves,  which  arise  from  the  pole  at  A,.  lying  upstream 
of  the  source  with  the  resultant  propagation  of  energy  toward  the  source.  Haskind  [4]  (as  does 
Bougis  [  1 ))  uses  a  fictitious  viscosity  n,  which  is  ultimately  set  equal  to  zero,  which  serves  to  give  the 
proper  placement  of  the  contours  around  the  poles.  He  shows  that  the  contour  should  always  pass 
above  k2  but  the  contour  around  k ,  depends  on  8.  In  particular,  for  cos  8  >  0,  i.e..  8  <  n/ 2.  the 
contour  should  pass  below  the  pole  k,.  With  this  change,  the  wave  patterns  obtained  by  Sretenskii 
would  be  similar  to  those  of  Becker  111),  who  shows  two  sets  of  wave  patterns  for  "  $  1/4. 


Bougis,  who  carries  out  a  systematic  and  detailed  implementation  of  the  above  procedure,  shows 
that  there  are  actually  three  regions:  0<  v  <  1/4,  1/4  <  v  <  1/^2,  and  I/v^  <  v.  These  pat¬ 
terns  are  reproduced  in  Fig.  3.  This  figure  shows  that  for  v  <  1/4  there  are  five  waves,  a  set  of 
ring  waves  and  two  complete  sets  of  Kelvin  like  waves,  each  composed  of  transverse  and  divergent 
waves.  For  v  >  1/4,  there  are  four  waves  with  a  region  in  front  of  the  source  which  is  now  distur¬ 
bance  free.  The  ring  waves  and  the  inner  transverse  and  divergent  waves  remain,  but  the  outer  set  of 
transverse  waves  has  disappeared^  leaving  only  the  divergent  waves.  The  only  difference  between  the 
cases  1/4  <  v  <  1  /^2  and  1  /'^2  <  v  is  that  in  the  former  case,  the  forward  part  of  the  ring  waves, 
near  the  wedge  lines  defining  the  boundaries  of  the  disturbance  free  region,  has  0-values  which 
correspond  to  forward  propagating  waves.  This  fact  had  previously  been  noted  by  Eggers  [12], 


3.2  Method  of  Stationary  Phase  Applied  to  Convolution  Time  Domain  Approach 


Brard  [6,7]  also  uses  the  method  of  stationary  phase.  However,  instead  of  the  single  application 
noted  above,  this  method  must  be  applied  three  times  to  the  triple  integral  time  domain  Green’s  func¬ 
tion  given  by  Eq.  (2.29):  first  with  respect  to  0,  then  with  respect  to  k ,  and  finally  with  respect  to  t. 
This  last  application  leads  to  a  solution  of 


d^l(T)ldT  =  0 


(3.8a) 


with  'k(r)  given  by 


*(r)  = 


’At  ~  T) 
4  Rr 


T  03  T 


(3.8b) 


where  Rr  is  defined  in  Fig.  4.  This  figure  shows  a  sketch  of  a  triangle  whose  vertices  are  the  field 
point  m ,  and  sT  and  s,,  the  projections  of  the  positions  of  the  source  at  times  t  and  t  onto  the  undis¬ 
turbed  free  surface.  The  lengths  of  the  sides  are  Rr,  R,  and  U(t  -  r),  and  the  base  angles  of  the  tri¬ 
angle  may  be  interpreted  as  being  y  and  0,  respectively  the  field  point  angle  and  wave  propagation 
direction,  as  above.  Upon  noting  that 


=  -  U  cos  0 


dV(T)/dT  =  0  leads  to 


cos  0  = 


U(t  -  T) 


U\t  -  t)2 


(3.10) 


Using  the  relation 


(3.11) 


U(t  -  t)  sin  (0  +  y) 

converts  Eq.  (3.10)  to  an  equation  between  y  and  0,  similar  to  Eq.  (3.7).  Physically,  Eq.  (3.10) 
corresponds  to  the  finding  of  the  positions  of  the  source  at  times  t,  st ,  which  give  waves  which  reach 
field  point  m  at  current  time  t . 


3.3  Geometrical  Construction  Technique 


Becker  considers  the  case  of  a  fixed  pulsating  source  in  the  presence  of  a  uniform  stream  U  in 
the  -x  direction.  He  starts  by  noting  that  for  an  observer  fixed  in  space,  all  wave  processes  occur  at 


1 


the  source  oscillation  frequency  /.  Also,  relative  to  a  fluid  at  rest,  the  wave  crests  propagate  in  the 
direction  0  at  phase  velocity  V ,  which  is  a  function  of  wavelength  X.  In  the  presence  of  the  current 
(/,  which  retards  the  wave  motion  by  the  velocity  U  cos  0,  the  resultant  velocity  is  therefore 
V  -  U  cos  0.  The  frequency  /  of  waves  propagating  past  a  fixed  point  is  then  related  to  this  vel¬ 
ocity  by 


m  -  t/ c <*8  _  f  ,3., 2) 

X 

Since  /  must  always  be  positive,  the  signs  in  the  numerator  must  be  switched  if  E(X)  <  U  cos  0. 
With  this  understanding,  Eq.  (3.12)  may  be  solved  for  E(X),  as  follows 

V  cos  0  ±  f  X  =  E(X)  (3.13) 

where  the  upper  or  lower  sign  is  taken  according  to  whether  U(\)  is  greater  or  less  than  V  cos  0. 
The  intersections  of  the  two  straight  lines 

/ 1  (X)  =  U  cos  0  +  f  X  (3.14a) 


f2  (X)  =  U  cos  0  -  /  X  (3.14b) 

with  the  K(X)  curve  gives  the  wavelengths  as  a  function  of  6,  X(0). 

In  order  to  obtain  a  6  -  y  relation,  similar  to  the  preceding  Eqs.  (3.7)  and  (3.10),  Becker  uses 
the  following  two  properties  of  the  n  th  curve  of  constant  phase 

4>  (jc  ,y;0)  as  nX(0)  -  xcos  0  -  y  sin  0  =  0  (3.15) 


d<t> 

00 


+  x  sin  0  -  v  cos  0=0. 
d6 


(3.16) 


An  expression  for  d\/dO  is  obtained  by  using  Eq.  (3. 13).  It  turns  out  that  the  same  0  -  7  relation  is 
obtained  by  using  either  the  upper  or  lower  sign.  Differentiating  Eq.  (3.13)  (with  the  upper  sign) 
with  respect  to  0  results  in 


-U  sin  0 


dV  d\  d\ 
d\  dO  J  dd  ' 


(3.17) 


Using  the  following  relationship  between  the  phase  velocity  V  and  the  group  velocity  U  to  obtain  an 
expression  for  dV !d\ 


converts  Eq.  (3.17)  to 


K(X)  - 


d\  _  U  sin  0 
d$  V  -  V 


(3.18) 


(3.19) 


Using  Eqs.  (3.12)  and  (3.15)  to  eliminate  /  and  X  from  the  right  hand  side  of  the  above  equation 
results  in 


dd  n  V„  -  U  cos  8 


( x  cos  0  4 -  y  sin  8). 


(3.20) 


Substituting  this  equation  into  Eq.  (3.16)  and  expressing  x  and  y  in  terms  of  the  polar  coordinates  R 
and  y  results  in  the  desired  y  -  8  relationship 


tan  (y  -  8)  = 


U  sin  8 
Vg  -  U  cos  8 


(3.21) 


which  may  also  be  transformed  to  the  more  direct  form 


tan  7  = 


Vg  sin  6 
Vg  cos  8  -  U 


(3.22) 


For  the  deep  water  case  considered  in  this  report,  Vg  is  given  by 

V  =iy=iJ^ 

*  2  2^  It' 


(3.23) 


Becker  then  considers  a  velocity  triangle  whose  sides  are  0  directed  along  the  x  axis,  Vg  at 
angle  6,  and  the  resultant  energy  transport  VE  given  by 


VE  =  U  +  Vg. 

He  shows  that  the  angle  which  VE  makes  with  the  x-axis,  y* ,  is  given  by 


(3.24) 


tan(y*  -  8)  = 


U  sin  8 
Vg  -  U  cos  8 


i.e.,  comparing  with  Eq.  (3.21) 


y*  =  y_ 


(3.25) 


(3.26) 


That  is,  at  a  given  point  {x  ,y)  in  space,  a  wave  normal  is  directed  at  angle  8,  but  the  resultant  energy 
is  transported  at  the  polar  angle  7  of  the  point. 

Eggers  [12]  starts  with  the  same  frequency  condition  as  Becker,  Eq.  (3.12).  However,  he  car¬ 
ries  out  his  analysis  in  terms  of  k  instead  of  X,  and  explicitly  introduces  the  deep  water  relation  (Eq. 
(3.23)  expressed  in  terms  of  k  =  2t/\)  from  the  start  instead  of  at  the  end.  As  a  result,  his  analysis 
is  of  a  more  mathematical  nature  and  to  a  large  extent  resembles  the  method  of  stationary  phase 
described  in  section  (3.1)  of  this  chapter. 

4.  ESTIMATES  OF  FAR  FIELD  WAVE  AMPLITUDE 

Only  Grekas  [2,3]  and  Haskind  [4]  give  explicit  formulas  for  the  amplitudes  of  the  far  field 
waves.  As  shown  below,  the  formulations  are  for  small  values  of  the  parameter  v  =  Uw/g  or  the 
velocity  U . 


4.1  Grekas’s  Approximate  Formula 

Since  Reference  [3]  is  basically  a  summary  of  Part  1  of  the  doctoral  thesis  [2],  the  description 
here  summarizes  the  derivation  given  in  Chapter  4  of  [2].  The  analysis  is  carried  out  in  a  fixed  coor- 


>1  M  Ul  <.l  «,t  t^i  «. 


dinate  system  for  a  pulsating  source  fixed  in  the  presence  of  a  horizontal  stream.  Grekas  actually 
considers  the  stream  to  have  components  in  both  horizontal  directions  x  and  y.  To  conform  to  the 
convention  used  in  this  report,  with  no  actual  loss  of  generality,  his  analysis  will  be  given  for  flow  U 
in  only  the  +x  direction.  The  crucial  simplification  of  his  approach  is  to  consider  U  small  such  that 
U2  can  be  neglected  everywhere.  In  the  case  of  the  free  surface  condition,  Eq.  (2.13)  is  simplified  to 

Taking  the  strength  of  the  pulsating  source  to  be  given  by 

Q  =  Q*  cos  wt  (4.2) 

G  is  then  decomposed  into  the  usual  sum  of  a  Rankine  source  and  a  function  which  is  harmonic 
everywhere  in  the  lower  half  space 

G  -  G,  +  G2  -  - — -  +  G,  (4.3) 

Putting  G2  in  the  usual  Fourier  transform  form,  Eq.  (2.6),  and  again  noting  the  key  identity  (2.7),  the 
use  of  the  free  surface  condition  (4.1)  leads  to  G  in  the  form 


-Q*  cos  ot 
4ir  r 


,*(z  +  z'  +w) 


- Re  1  e-‘«  \  de  l  - f - BL2* -  (4.4) 

4X2  [_  -(w2  -  2kuU  cos  6  -  gk)  J 

where  w  is  defined  in  Eq.  (2.7b).  The  neglect  of  U2  in  the  free  surface  condition  simplifies  the 
denominator  of  the  integrand  to  gk  +  2k uU  cos  6  -  to2,  which  has  a  simple  pole  in  k ,  whereas  the 
more  general  free  surface  condition  (2.2)  leads  to  the  more  complex  denominator  in  Eq.  (2.12), 
gk  -  («  T  Uk  cos  0)2,  which  has  a  double  pole  in  k . 

The  equation  for  the  “meridian”  lines  (or  wave  crest  contours)  of  the  deformed  free  surface  is 
taken  to  be 

h  =  —  it  =  h'  +  h*  <45) 

g  ot 

where  h,  is  a  near  field  solution  which  goes  to  zero  as  l/Rn ,  n  >  1,  and  h2  is  the  far  field  solution. 
It  should  be  emphasized  h  is  not  the  total  wave  elevation  since  this  would  include  the  convective  term 
(U/g)d<i>/d x.  This  term  can,  however,  be  easily  obtained  from  the  final  form  for  h2 ,  given  below  in 
Eq.  (4.9).  The  far  field  solution  is  put  in  the  usual  form 


h7  =  - 


Re  {  e  ~iwl  f(6)  e,R*t9)  d6 

2  x 


ijV 

>»w 


where 


g(d)  =  k0  cos  (d  -  0)  = 


or  cos  (6  -  0) 
g  +  2o)U  cos  6 


m  = 


g  +  2co£/  cos  6 


0  =  tan  1  y  /x  is  the  spatial  angle. 


The  quantity  it  0  is  the  wavenumber  in  terms  of  the  intrinsic  frequency  w0 

wo  (w  -  V  cos  6  k 0)2 

k  -  —  =  - 


fa)2  -  2u  U  cos  6  k0 


(4.7a) 


(4.7b) 


where  Eq.  (4.7b)  results  by  again  neglecting  the  U2  term.  Solving  Eq.  (4.7b)  gives  ik0  in  terms  of  u 
as  follows 


g  +  2  wt/  cos  0 

Applying  the  method  of  stationary  phase  to  Eq.  (4.6),  the  following  final  form  for  h2  is  obtained 


/<w 


2l r  il/koWo)/?  cos  (@o-0)-*/4-u>r) 

l*"(*o)l  * 


where 


,,,« .  =  zl cos  {6°  ~ g)  0,2 

0  (g  +  2w{/  cos  60)2 

and  0O  is  the  root  of  the  6-0  relationship  given  by 

g'0 0)  =  (P  ~  #o)  +  2lj)U  sin  0-0. 


(4. 10a) 


(4.10b) 


Grekas's  expression  for  g”(0o)  erroneously  has  the  denominator  raised  to  a  power  of  3  instead  of  2. 
Figure  5  shows  that  there  is  only  a  set  of  ring  waves,  as  opposed  to  the  ring  waves  plus  two  sets  of 
ship  transverse  and  divergent  waves  shown  in  Fig.  3  for  v  <  1/4. 

4.2  Haskind’s  Approximate  Formula 

Haskind  (4)  uses  a  moving  coordinate  system,  for  which  the  free  surface  condition  is  given  by 
Eq.  (2.2).  The  Fourier  transform  technique  is  also  used  to  obtain  G2  in  the  expression 


Ge‘*  =  -  -  —  +  G,  eiul. 

r  r\ 


(4.11) 
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A  fictitious  viscosity  n  is  included  in  his  formulation.  This  viscosity,  although  ultimately  set  equal  to 
zero,  serves  to  indicate  the  proper  way  of  surrounding  the  poles  in  the  integration  with  respect  to  k . 
Basically,  the  semi-circles  surrounding  the  poles  must  be  chosen  such  that  their  imaginary  parts  must 
have  the  same  sign  as  the  imaginary  parts  of  the  poles  k\,  k'2  when  ^  *0 


*1,2  =  *1,2  “  if  W 

where  *i,*2  are  given  in  Eq.  (3.5b).  He  shows  that 


sign  Im  k  j  =  -  sign  cos  8 


Uu 


sign  Im  =  -  sign  v  —  —  sign - . 

8 


integration  for  the  integral  with  respect  to  k 

G2  =  F, (x,y,z)  +  F2  ( x,y,z )  +  F3  (x  ,y  ,z) 


(4.14) 


where 


1  ,*/2  .<*>  k  ,,(*(* +*'+'w>l 

F,  = - -  (  f  - - -  - - dk  d$ 

r  * -*/2  (jt  —  fc,)  v  1  +  4  v  cos  8 


J_  f ,i/2  k  e^+t'-iw)] 

T  ‘  ~T/2  W,)  (k  -  *,)  Vi  -  4  V  cos  8 


dk  d8 


F  ~  —  r+T_V 

2  T  J-*  +  »o  %.,) 


falkU+z'+iw)  1 


K  »+90 

-  -4  L. 


(fc  -  *2)  ^1  +  4i>  cos  8 

o  ^|*(Z+Z  +lV)] 


dk  dO 


tv2  J*-«o  0  (^  _  £,)  (k  -  k2)  cos2  8 


dk  d8 


*1,2  (0)  -  *1.2  (0  +  T) 

'  signifies  that  the  integration  is  to  exclude  the  interval  (-0o,  +60) 
AT  =  «2/$ 
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(4.12) 
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The  contour  L  j  connects  the  points  k  =  0_  and  k  =  oo  by  passing  above  k  =  it , ,  the  contour  L ,  con¬ 
nects  these  points  by  passing  below  k  =  kt,  and  die  contour  L2  passes  above  k  =  k2. 


For  values  of  |  v  | 


Haskind  gives  the  following  asymptotic  far  field  formula  for  G2  in 


terms  of  the  polar  coordinates  R ,  0 


G2  =  2(1  ±  2  v  cos  0) 


,1*0  +2v  cos  0)1;  +  ;'  +  HR  -x‘  cos  0-v’  sm  0)]  +  /»/4| 


(4.15) 


where  the  upper  sign  is  to  be  used  for  positive  values  of  v  and  the  lower  sign  for  negative  values  of 
v.  Unfortunately,  Haskind  does  not  provide  any  details  of  the  derivation.  However,  similar  to  the 
formulation  by  Grekas,  Eq.  (4.9),  the  above  equation  gives  rise  to  a  single  system  of  ring  waves 
which  propagate  in  all  directions  from  the  source.  This  again  suggests  that  the  free  surface  condition 
(2.2)  gives  rise  to  a  simple  pole  in  k.  Since  the  term  U2  d2/dx2  gives  rise  to  a  square  term  in  k,  it 
must  have  been  neglected  in  the  derivation. 
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